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In spite of increasing numbers of countries having established renewable energy development 
mechanisms for carbon dioxide (CO2) emissions reduction, the COz emissions problem continues to 
worsen along with the growth of the world economy. This leads us to examine the threshold effect of the 
proportion of renewable energy supply for COz emissions reduction by means of the panel threshold 


Keywords: regression model (PTR). Economic growth and the price of energy are also both taken into account in the 

Renewable energy model in measuring the specific influence that each of them has on CO, emissions. The empirical panel 

=e oa data encompass all 30 member countries of the OECD and cover a period of about a decade in length from 
ea 


1996 to 2005. Our empirical results provide clear evidence of the existence of a single threshold effect 
that may be divided into lower and higher regimes. Based on the specific estimates of the slope 
coefficients in each regime distinguished, we find that a renewable energy supply accounting for at least 
8.3889% of total energy supply would mean that CO, emissions would start to be mitigated. Furthermore, 
real GDP and the CPI of energy are significantly and positively and insignificantly and negatively 
correlated with CO, emissions, respectively. These findings lead us to conclude that the authorities ought 
to enhance the proportion of renewable energy supply to more than 8.3889% of all energy supplied, 
which might help resolve the dilemma between economic growth and CO, emissions. Realizing the 
effects of CO emissions reduction via energy price reforms or the levying of a carbon tax levy may, 
however, still remain a puzzle. 
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1. Introduction 


The Kyoto Protocol, which came into force on February 16, 
2005, places obligations on all signatories to ensure that 
greenhouse gas emissions in 2012 are no greater than the total 
of such emissions in 1990. Many of the countries that seek to 
achieve this goal focus on three main strategies that are 
available to them for decreasing the amount of carbon dioxide 
(CO2), namely, the greater efficiencies in electricity generation, 
as well as the recycling, capturing, utilization, disposal/storage, 
and use of renewable and nuclear sources of energy [1]. These 
strategies also suggest that fossil energy use be diminished 
based on well-established substitution by renewable energies. 
Nevertheless, most renewable sources are weaker economically 
than fossil fuel sources. This is a reason why many countries 
provide incentives (e.g., feed-in-tariffs, investment subsidies 
and tax credits, etc.) to enhance the proportion of their 
renewable energy supply. 

Renewable energy gives rise to fewer CO, emissions than fossil 
energy such as biomass, and even zero emissions in the case of 
wind power, solar photovoltaic, etc. A substantial reduction in CO3 
emissions can be achieved through a large-scale integration of 
renewable energy supply into the energy system. Previous studies 
simulate optimal scenarios of the joint linear or nonlinear effects 
of renewable energy supply on CO2 emissions reduction. The 
results describe the positive relationship between renewable 
energy supply and CO2 emissions reduction after renewable 
sources are included in a portfolio for the energy system [2-9]. Itis 
found that renewable energy has developed rapidly with a 
consensus being reached in CO; emissions reduction and with the 
high volatility in crude oil prices in recent years. The global 
temperature continues to increase year by year as usual. This 
phenomenon has led to increased interest in reestimating the 
specific effects both below and above the optimal threshold 
parameter of the proportion of renewable energy supply for CO2 
emissions reduction. Furthermore, we investigate the individual 
relationships between CO2 emissions, economic growth and the 
energy price. These latter two are proxied by real gross domestic 
product (GDP) and the consumer price index (CPI) of energy, 
respectively. 

In the current academic literature, there is evidence that with 
more economic growth, energy consumption has increased. The 
significant positive relationship between economic growth and 
CO- emissions has been clearly observed [10-15]. One strand of 
the literature has investigated the existence of a causal relation- 
ship between the greenhouse gas emissions and economic growth 
indices using multiple tests on climate change such as the famous 
Environmental Kuznets Curve which summarizes the relationship 
as taking the form of an inverted U-shape [16-18]. Another line of 
investigation in terms of the causality issue of economic growth 
and CO, emissions is based on time series or panel data 
econometric techniques, such as unit root tests, cointegration 
and the related error correction model. The nature of the causality 
between per capita GDP and per capita CO2 emissions, short-run 
dynamic comovements and long-run cointegration has been 
explored [19]. These conclusions in the literature would appear to 
suggest that the authorities ought to exercise caution due to the 
dilemma that exists between economic development and 
environmental protection in terms of the long-run balanced 
relationship. 

The debate on CO» emissions reduction in light of the steep 
increases in energy prices is never-ending. Intuitively, the energy 
consumption will decline as energy prices are raised, thereby 
resulting in a decrease in COz emissions. High energy prices lead to 
a significant reduction in the total primary energy supply and this 


is combined with structural changes in primary energy supply in 
view of the smaller quantities of oil and natural gas through a 
noticeable increase in renewable energy sources [20]. Energy price 
reform has thus become a crucial element in the promotion of 
energy conservation, a reduction in CO3 emissions and the 
substitution of fossil fuels with renewable energy. Deregulating 
energy prices and imposing different levels of taxation on fossil 
fuels could decrease CO emissions without considerably suppres- 
sing the growth of the economy. However, these policies cannot 
cope with the inevitable increased use of coal in the long run [21]. 
The literature thus induces us to hypothesize that an uncertain 
negative relationship exists between energy prices and CO, 
emissions. 

We examine the panel threshold regression model (hereafter, 
PTR) that was previous described, whose estimation methodol- 
ogy is developed for non-dynamic panels with individual specific 
fixed effects. Least squares estimation of the threshold and 
regression slopes is proposed using fixed-effects transforma- 
tions. A non-standard asymptotic theory of inference and a 
bootstrap method are also imposed which allows for the 
construction of confidence intervals to assess the statistical 
significance of the threshold effect and the testing of hypotheses 
[22]. The PTR model can then estimate the optimal threshold 
parameter which is divided into lower and higher regimes. There 
is a specific estimation of the slope coefficient in each regime. 
Our empirical panel data set is formed across all 30 member 
countries of the Organisation for Economic Co-operation and 
Development (OECD)! and the sample period spans a decade 
from 1996 to 2005. 

The most important contribution of this paper is that we are 
able to determine the optimal proportion of renewable energy 
supply for CO2 emissions reduction. Our empirical results show 
that only the single threshold effect is explored. The estimated 
specific slope coefficients for the lower and higher regimes have a 
completely opposite significance to those for CO2 emissions. We 
might interpret this robust evidence as indicating that the 
proportion of renewable energy supply is at least more than 
8.3889% (of the estimated optimal threshold parameter) that 
would initially mitigate CO emissions. 

Other empirical results also show that real GDP is both 
significantly and positively correlated with CO, emissions. This 
relationship implies that economic growth will increase the 
environmental impact as global temperatures abnormally rise. 
An analogous financial instrument is a carbon tax which has a 
similar effect in terms of raising energy prices for mitigated CO, 
emissions [23-25]. There is an insignificant negative relationship 
between the CPI of energy and CO2 emissions. Hence realizing the 
effects of a reduction in CO, emissions via energy price reform or a 
carbon tax levy may still be a puzzle. To sum up, our empirical 
findings lead us to conclude that the authorities ought to enhance 
the proportion of renewable energy supply by more than 8.3889%, 
if the dilemma between economic growth and CO; emissions is to 
be resolved. 

The remainder of this paper is organized as follows. Section 2 
describes the PTR model estimation methodology. Section 3 
describes how the panel data are formed. Section 4 presents the 
results of the estimation which are robust with respect to the 
threshold effect. The final section summarizes our empirical 
findings and enables us to draw conclusions. 


1 member countries of the OECD are Australia, Austria, Belgium, Canada, Czech 


Republic, Denmark, Finland, France, Germany, Greece, Hungary, Iceland, Ireland, 
Italy, Japan, Korea, Luxembourg, Mexico, the Netherlands, New Zealand, Norway, 
Poland, Portugal, Slovak Republic, Spain, Sweden, Switzerland, Turkey, United 
Kingdom and United States. 
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2. The methodology 
2.1. Single threshold estimation 


The empirical data is from a balanced panel. Our empirical 
objective is to examine whether there exists a renewable energy 
supply threshold effect for a reduction in CO, emissions. In order to 
present the notation, let the dependent variable CO2; be the 
annual growth in terms of the percentage of CO, emissions, the 
independent variable GDP;, be the annual growth in terms of 
the percentage of real GDP, CPI, be the annual growth in terms of 
the percentage of the CPI of energy, and RES; be the annual 
contribution of renewables to energy supply. The single threshold 
estimation of the PTR model as follows: 


CO2i¢ = pi + By GDPit + By CPlit + B3 RESit I(RESit 
<y)+ By RESjt I(RESjt > V) + eit (1) 


where intercept term u; is implied the individual specific mean. 
The subscript i indexes the individual, the subscript t indexes times 
and I( ) is the indicator function. GDP;,, CPI; and RES; are used to 
measure the specific influences on CO2;, simultaneously. The other 
RES; in the indicator functions I( ) is the threshold variable which is 
used to estimate the optimal threshold parameters 7. There are 
diverse into lower or higher regimes for which the specific 
estimates of the slope coefficients are 6, and 6, for each regime 
distinguished. The error term e;; is assumed to be independent and 
identically distributed (iid) with mean zero and finite variance o°. 


2.2. Estimation 


An alternative compact representation of (1) is to set RES;(yv) = 


RESjz I(RESit < y) rg I 
Eo IRES: > y) and = (3 4) so that (1) equals 


CO2; = pi + By GDPi + By CPli + B’ RESit(Y) + eit (2) 


One traditional method to estimate (2) that to remove 
individual effect u; that taking averages (1) over the time index 
t produces 


CO2; = u; + 61 GDP; + By CPI; + B’ RESi(y) + & (3) 


Taking the difference between (2) and (3) for de-mean that 
yields 


CO2}, = B GDP} + Pa CPI + p' RES}(Y) + ei (4) 


For any given y, the slope coefficient 6 can be estimated by 
ordinary least squares (OLS). The sum of squared error is calculated 
by 


Si(y) = Ë (Y)/E(y) (5) 


Estimation of y by least squares is easiest to achieve by 
minimization of the concentrated sum of squared errors S4(y) 
[26,27]. Hence the least squares estimators of y is 


y = argmin S41 (y) (6) 
Y 


Once Ŷ is obtained, the slope coefficient estimate is B = Bip). 
The residual vector is 6° = ê (f) and residual variance 


6? = E&E /n(T — 1) = S1(P)/n(T — 1) (7) 


where n indexes the number of sample, T indexes the period of 
sample. 


2.3. Testing for single threshold 


There is important to determine whether the threshold effect is 
statistically significant. The null hypotheses of no threshold effect 
HO: 63 = 64 in (1). Under this null hypothesis of no threshold, the 
model (1) is reduce as 


CO2it = jt Bi GDPit t Bs CPlit t Bs; RESit + Cit (8) 


The regression parameter £3 is estimated by OLS, yielding 
coefficient £3, residuals é;, and sum of squared errors Sp = & &. 
The likelihood ratio test of Hg is base on 


So — Sı (Ŷ 
0 FP (9) 


Fi(y)= 

The asymptotic distribution of Fi(y) is non-standard, and 
strictly dominates the X? distribution. A bootstrap procedure 
attains the first-order asymptotic distribution which the null of no 
threshold effect is rejected if F,;(y) is higher than the desired critical 
value [28]. 

When there is a threshold effect (61 4 62), that y is consistent 
for yo (the true value y). The best way to form confidence intervals 
for y is to form the no-rejection region using the likelihood ratio 
statistic for tests on y [27]. The null hypotheses is HO: y= yo, the 
likelihood ratio test is to reject for large values of LR;(yo) where 


The theorem LR, y — dé, as n — oo, where & is a random variable 
with distribution function 


P(E<x) = (1- exp(-2)) (11) 


The asymptotic distribution of the likelihood ratio statistic is 
non-standard yet free of nuisance parameters. The technical 
assumptions include the rather unusual condition that 64 — 63 — 0 
as n— co. The asymptotic distribution may be used to valid 
asymptotic confidence intervals. The distribution function (11) has 
the inverse 


c(a@) = —2log(1 — v1 —a@) (12) 


from which calculate critical values of rejects at the asymptotic 
level. For example, the 10%, 5% and 1% critical value are 5.94, 7.35 
and 10.59, respectively. 


2.4. Double threshold estimation 


Extend model (1) has been double threshold model takes the 
form 


CO2i¢ = pi + By GDPi + By CPlit + B3 RESit IRES < 71) 
+ Bi RESit (y; < RESit < V2) + Bs RESit (V2 < RESit) 
+ Cit (13) 


Thus for given (y1, y2) the concentrated sum of squared errors 
S(yı, y2) is straightforward to calculate as the single threshold 
model. Let S;(7) be the single threshold sum of squared errors as 
defined in (5) and let ~, be the threshold estimate which 
minimizes S;(y). The 7, will be consistent for either yı or yo. 
Fixing the 7,, then criterion S}(y) is 


[Siy if Pi<y 
S5 5 1 v2 5 1 v2 14 
2(Yo) ees if y.<Y (14) 


1672 C.-L. Chiu, T.-H. Chang/ Renewable and Sustainable Energy Reviews 13 (2009) 1669-1674 


and threshold estimate 7} is 


7, = argmin S} (72) (15) 
Y2 
Fixing the asymptotic efficiency of f, define the refinement 
criterion 


r = S(1, P2) if 1<% 
Son) = {Sor if hey, i 


and the refinement estimate 
ĵi = argmin S{ (y1) (17) 
Yı 


The minimizing sum of squared errors S}(f,) with variance 
estimate ô? = S}(f3)/n(T — 1). Thus an approximate likelihood 
ratio test of double thresholds can be base on the statistic 

S1 (P1) — SP 
Fa(y) = MPa) — 2(V2) (18) 
ô 
and construction of confidence intervals in the same way as 
previous Section 2.3 for the individual threshold parameter as 


S) - S594) 


LR5 (y) = and L(y) ==) gy 


6 

The computations of the least squares estimate of the threshold 
parameter y, yı and y> involve the minimization problem (6), (15) 
and (17). Since the sum of squared error function S;(y), SY (y4) and 
S¥(y2) depends only through the indicator functions are step 
function with at most nT steps, with the steps occurring at distinct 
values of the observed threshold variable. Thus the optimization 
search describe above may be numerically intensive. One approach 
may be taken. Sort the distinct values of the observations on the 
threshold variable. Instead of searching over all values of threshold 
variable the search may be limited to specific quantiles. This greatly 
reduces the number of regressions performed in the search. For the 
empirical work, we use the grid {0.25%, 0.50%, 0.75%, 1.00%,..., 
99.25%, 99.50%, 99075%, 100%} which contains 400 quantiles. 


3. The data 


The panel data set used in this study encompasses all 30 
member countries of the OECD and the sample period extends over 
a decade from 1996 to 2005. The data are obtained from the OECD 
statistical database, and all variables are defined at the end of the 
year. One of the variables is the contribution of renewables to 
energy supply. The other variables are the specific growth in the 
percentage of CO, emissions from energy use, real GDP and the CPI 
of energy. It should be noted that we calculate the annual growth in 
percentage of CO2 emissions, real GDP and the CPI of energy based 
on the logarithmic difference. The total number of observations for 
the cross-sectional and time series data is 300 for each variable. 

Table 1 exhibits the summary statistics for each of the four 
variables. The sample mean of the annual growth in the percentage 
of CO; emissions from energy use is found to be positive and to be 
1.1625%. This statistic indicates that the global greenhouse effect 


Table 1 
Summary statistics. 


problem is continually being exacerbated by CO2 emissions. On 
average, the annual contribution of renewables to energy supply is 
12.1777%, and the standard deviation is 14.8662 or higher, which 
means that the development of renewable energy differs quite 
significantly among the member countries of the OECD. The 
maximum contribution of 73.8292% is found in Iceland in 2005, 
which means that this country has the highest renewable energy 
supply among the member countries of the OECD, with its average 
contribution exceeding 70% from 1996 to 2005 in each year. In fact, 
large-scale renewable sources of hydropower and geothermal 
power have been exploited in Iceland. In contrast, the minimum 
annual contribution of renewables to energy supply of 0.8095% is 
found in Korea in 1996. Only slight growth, as evidenced by the 
contribution of renewables to energy supply reaching 0.9615%, is 
found to have occurred a decade later. 

Table 1 also reports that the sample mean of the annual growth 
in percentage terms of real GDP is positive and amounts to 
3.1979%, having a standard deviation of 2.3053. Furthermore, the 
sample mean of the annual growth in terms of the percentage of 
the CPI of energy is positive and amounts to 6.3876%, with a 
standard deviation of 10.1142. Both statistics imply that the rapid 
rise in energy prices is accompanied by higher volatility. The 
maximum such price rise is 72.6484% in Turkey in 1996. This is the 
highest energy price growth recorded among the member 
countries of the OECD, for which the average fluctuation exceeded 
40% in each year from 1996 to 2005. The Jarque-Bera statistical 
tests indicate that all variables significantly reject the normal 
distribution hypothesis at the 1% level, which indicates that the 
distribution of each time series is both leptokurtic and fat-tailed. 


4. Empirical results 
4.1. Threshold effect estimation and test 


To determine the optimal threshold parameters from Eq. (13), 
we use least squares estimation, and allow for zero, one and two 
thresholds. The test statistics of F,(y) and F2(y) along with the 
bootstrap critical values are reported in Table 2. The test for a single 
threshold F,(y) is 22.5729. There is a greater significance level than 
the 95% critical value which is 22.4047. The test for a double 
threshold F,(y) is 18.4303. There is a significance level that is less 
than the 90% critical value that is 18.7472. These results provide 
strong evidence of only a single threshold effect being uncovered in 
the regression. 

The point estimates of the threshold asymptotic 95% confidence 
intervals are also reported in Table 2. The estimation of 7, yields a 
value of 8.3889 in the empirical distribution of the threshold 
variable. Thus the proportion of renewable energy supply has been 
separated into lower and higher regimes. The asymptotic 
confidence intervals for the threshold 7, are very tight and lie 
within the range of 8.3889 to 8.3889, indicating that there is little 
uncertainty regarding the nature of this breakdown. The estima- 
tion of 7, yields a higher value of 11.3122, for which the 
asymptotic confidence intervals are very wide and quite uncertain 
in that they range from 10.6473 to 23.3276. The findings are robust 
in that the lower regime I(RES; < ~,) includes 182 (60.67%) 


Variable 


Mean (%) (standard deviation) 


Minimum (%) Maximum (%) Jarque-Bera (p-value) 


Annual growth in percentage of CO emissions from energy use 
Annual contribution of renewables to energy supply 

Annual growth in percentage of Real GDP 

Annual growth in percentage of CPI of energy 


1.1625 (4.6680) 
12.1777 (14.8662) 
3.1979 (2.3053) 
6.3876 (10.1142) 


—21.7716 20.3322 126.2835 (0.0000) 
0.8095 73.8292 749.0096 (0.0000) 
—7.4954 11.6801 132.0706 (0.0000) 
—11.0988 72.6484 3341.9191 (0.0000) 
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Table 2 
Threshold estimates. 


Threshold variable 


Parameter (95% confidence interval) 


Fy) (10%, 5%, 1% critical values) 


Îi 8.3889" (8.3889, 8.3889) 


Po 11.3122 (10.6473, 23.3276) 


22.5729 (19.7475, 22.4047, 26.1675) 
18.4303 (18.7472, 20.1217, 22.7566) 


Note: The critical values are simulated by bootstrap method. 
” Significance at the 5% level. 


observations and the relatively higher regime I(RES; > 7) includes 
118 (39.33%) observations. These results reveal that the single 
threshold model conveys information that is only found in the 
regression and the optimal threshold parameter is 8.3889%. 

Figs. 1 and 2 depict the construction of the confidence intervals 
for the single and double threshold models, respectively. The 
horizontal axis denotes the estimation of the first and second 
threshold parameters and the vertical axis expresses the concen- 
trated individual likelihood ratio function LR (y) and LR}(y). The 
point estimates are the values of 7, and 7, where the likelihood 
ratios LR} (y) and LR} (y) intersect the zero horizontal axis, which is 
in the far left part of the graph. The 95% confidence intervals for yı 
and y2 can be found from LRY (y) and LR} (y) based on the values of y, 
and y2 for which the likelihood ratio lies beneath the dotted line at 
7.35. This critical value is calculated using formula (12). 


4.2. Regression estimation 


According to the specification in the preceding section, a 
single threshold effect is found to exist. Thus the remainder of 


25 T T T T T T T 


Likelihood Ratio 


0 10 20 30 40 50 60 70 80 
First Threshold Parameter 


Fig. 1. Confidence interval construction in single threshold model. 
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Fig. 2. Confidence interval construction in double threshold model. 


Eq. (13) may be reduced according to the class of model (1). 
Where the regression slope estimates contain conventional 
OLS standard errors and White-corrected standard errors 
which are reported in Table 3. Only the coefficient fz is 
insignificant, and 6, is 0.4630 at the 1% significant level. There is 
an obvious positive relationship between real GDP and CO2 
emissions as the findings indicate in the previous literature, 
where an expansion in real GDP results in increased CO2 
emissions. More fossil fuels are consumed as economic growth 
increases, as fossil fuels are more stable and also more efficient 
than renewable sources of energy. Moreover, fossil fuels also 
provide cheaper electricity and fuel that further encourage 
economic growth. The coefficient $2 is very small and negative 
as shown by its value of —0.0002, which indicates that an 
uncertain negative relationship exists between the price of 
energy and CO2 emissions. Consequently, as to whether rising 
energy prices via a demand/supply matching mechanism or a 
carbon tax levy would benefit a reduction in COz emissions is a 
question that still remains unanswered in spite of our empirical 
findings. 

Table 3 also indicates that the coefficients of primary interest 
are those to do with the annual contribution of renewables to 
energy supply. When the proportion of renewable energy supply 
in one of our regimes is equal to or lower than the optimal 
threshold parameter 7,, the specific estimation of the slope 
coefficient 63 is 1.1162 at the 5% significance level. An obvious 
positive relationship implies that, in spite of the increased 
proportion of renewable energy supply, there is no effect that 
will help mitigate CO» emissions which will continue to 
increase. In another of the regimes selected, the proportion of 
renewable energy supply is higher than the optimal threshold 
parameter 7,, with the specific estimation of the slope 
coefficient 84 being —0.9186 at the 1% significant level. There 
is an obvious negative relationship that leads to a reasonable 
explanation in that the increased economic growth leads to 
increased energy consumption. Thus a smaller proportion of 
renewable energy supply fails to generate a real substitutive 
effect of CO emissions for fossil energy. It is also interesting to 
note that only when the proportion of renewable energy supply 
is at least more than 8.3889% (the optimal threshold parameter) 
will CO, emissions start to be mitigated. If this is the case, we 
would like to suggest that the authorities increase the 
proportion of renewable energy supply to more than 8.3889%, 
which might help resolve the dilemma between economic 
growth and CO; emissions. 


Table 3 
Regression estimates: single threshold model. 


Variable Coefficient OLS standard error White standard error 
Bi 0.4630” 0.1256 0.1792 
po —0.0002 0.0366 0.0382 
Bs 1.1162" 0.4875 0.5638 
Ba 0.9186" 0.1893 0.2617 


” Significance at the 5% level. 
™ Significance at the 1% level. 
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5. Conclusions 


We are concerned with the phenomenon that renewable energy 
supplies may fail to resolve the greenhouse gas emissions problem 
even after renewable sources of energy are included in a portfolio 
of the energy system. Hence we recommend that the PTR model be 
used to examine the threshold effect of the proportion of total 
energy to be accounted for by renewable energy supply to give rise 
to a reduction in CO, emissions. Real GDP and the CPI of energy are 
also included in the model to measure their specific influences on 
CO2 emissions simultaneously. 

Our empirical approach involves using panel data which 
encompass all 30 member countries of the OECD over a sample 
period spanning a decade from 1996 to 2005. The empirical tests 
that examine only a single threshold effect are also explored. 
The threshold effect is divided into both lower and higher 
regimes. There is a complete reversal of the specific estimation 
of the slope coefficient between each distinguished regime. The 
real GDP and the CPI of energy are significantly and positively 
and insignificantly and negatively correlated with CO2 emis- 
sions, respectively. 

To sum up, our empirical findings lead to the conclusion that 
the demand for more fossil fuel consumption resulting from 
increased economic growth will speed up the generation of CO, 
emissions. Thus the authorities ought to increase the proportion of 
total energy sources accounted for by renewable energy supply to 
more than 8.3889%, which might help to solve the growing 
dilemma between economic growth and CO; emissions within the 
framework of the Kyoto Protocol. Realizing the effects of reducing 
CO2 emissions via an energy price reform or a carbon tax levy, 
however, still remains a puzzle. 
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